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Abstract. We study the large deviations of a simple noise-perturbed dynamical 
system having continuous sets of steady states, which mimick those found in some 
partial differential equations related, for example, to turbulence problems. The system 
is a two-dimensional nonlinear Langevin equation involving a dissipative, non-potential 
force, which has the essential effect of creating a line of stable fixed points (attracting 
line) touching a line of unstable fixed points (repelling line) . Using different analytical 
and numerical techniques, we show that the stationary distribution of this system 
satisfies, in the low-noise limit, a large deviation principle containing two competing 
terms: i) a "classical" but sub-dominant large deviation term, which can be derived 
from the Freidlin-Wentzell theory of large deviations by studying the fluctuation paths 
or instantons of the system near the attracting line, and ii) a dominant large deviation 
term, which does not follow from the Freidlin-Wentzell theory, as it is related to 
fluctuation paths of zero action, referred to as sub-instantons, emanating from the 
repelling line. We discuss the nature of these sub-instantons, and show how they arise 
from the connection between the attracting and repelling lines. We also discuss in a 
more general way how we expect these to arise in more general stochastic systems 
having connected sets of stable and unstable fixed points, and how they should 
determine the large deviation properties of these systems. 

PACS numbers: 02.50.-r, 05.10.Gg, 05.40.-a 
1. Introduction 

The dynamics of physical systems described by partial differential equations, such 
as those appearing in hydrodynamics, optics or quantum physics, are in many cases 
qualitatively similar to that of finite-dimensional dynamical systems, especially when 
an important dissipation mechanism is involved. When this is the case, instabilities, 
bifurcations, limit cycles and attractors are indeed often similar to their finite- 
dimensional counterparts (see, e.g., |l]). By contrast, when dissipative mechanism do 
not exist or are very small in partial differential equations, these may exhibit phenomena 
that have no counterparts in finite-dimensional systems. Examples of such phenomena 
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include non-dissipative relaxation and asymptotic stability |2j-|4], solitons, as well as the 
appearance of an infinite number of conserved quantities and infinite number of steady 
states |5}]8). 

In this work, we are interested in studying the rare events or large deviations |9j of 
systems possessing continuous sets of steady states which effectively act as attractors. 
Such sets of steady states are found in many dynamical equations having a non- 
canonical Hamiltonian structure, such as the 2D Euler equation, the Vlasov equation, 
magneto-hydrodynamic equations, and the shallow-water equations, to mention only 
a few examples; see [5||8] for more. In all of these examples, the consequence of the 
non-canonical structure is the existence of an infinite number of conserved quantities or 
Casimirs, which are responsible for the infinite (and continuous) set of steady states [5]|8] . 
Physically, these states are important because they can act as attractors, as has been 



found in experiments 10,11 and numerical simulations 12. Moreover, in some cases 



their attractive behavior can be explained theoretically [3||4j using arguments and 
methods based on statistical mechanics [8 ,13 15 . 



The problem that we address here is how the existence of a continuous set of 
steady states of a system influences its large deviation properties and how these 
properties compare with those of finite-dimensional systems. For the latter systems, 
two generic classes of systems have been considered from the large deviation point 
of view, namely: i) strongly dissipative systems with well-defined and disconnected 
attractors, as exemplified by the gradient dynamics of a Brownian particle in a potential 
(Kramers problem) [l6], and ii) weakly-perturbed Hamiltonian systems with added 
weak friction and noise. The large deviations of both of these classes are well known: 
they can be obtained using semi-classical (WKB or instanton) approximations of path 
integrals [9| [l7[]22] or, more rigorously, using the theory developed by Freidlin and 
Wentzell |23). 

Dynamical systems having an infinite number of steady states do not strictly fall in 
either of these generic classes. The fact that the dynamics of these systems is effectively 
irreversible and converge towards attractors that are asymptotically stable suggests that 
they are analogous to the class of strongly dissipative systems. However, because their 
steady states form one or several connected sets, the presence of a weak noise should lead 
these systems to diffuse over their attractors, as is the case for Hamiltonian systems. 
From this point of view, systems with infinite steady states share aspects of both strongly 
dissipative systems and weakly-perturbed Hamiltonian systems. 

To illustrate this point, we consider in this paper a simple model with two degrees of 
freedom, denoted by A and B, having a continuous set of steady states in the zero- noise 
limit. Our goal in the following is to obtain the stationary probability density P{A, B) 
of this model, which we refer to simply as the AB model, in the low-noise limit — )■ 0, 
where is the noise amplitude. Following the theory of large deviations, we expect 
this density to have the approximate form 

P(A5) ~e-'=''^(^'^), (1) 
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in the limit z/ — )■ 0, where Cu is a coefficient, called the speed, which diverges as z/ — )■ 0, 
and I{A, B) is a //-independent function, called the rate function or quasi-potential. 
Large deviation approximations of this form have been extensively studied, as mentioned 
above, for various noise-perturbed dynamical systems, including systems with single- 
point attractors as well as systems with attracting limit cycles (see, e.g., 17-22]). In 
nearly all the classical examples that we are aware of, the speed of the large deviation 
approximation is 1/z/. In the case, for example, of a Langevin dynamics describing 
the overdamped motion of a Brownian particle in a potential V{x) and in a fluid at 
temperature T, the speed is l/{kBT), with fc^ the Boltzmann constant, while the rate 
function is the potential V{x). The resulting large deviation form for the stationary 
density, 

P{x) ^ e-yi-)/ikBT)^ (2) 

is in agreement with the known Arrhenius factor describing transitions probability 
between stable or metastable states (Kramers theory), and implies, as physically 
expected, that P{x) concentrates in the low-noise limit (or low-temperature limit) on 
the equilibrium state minimizing the potential. 

The results that we obtain for the AB model show that P{A, B) concentrates in a 
similar way on the set of stable steady states of the model - in this case, a whole line of 
steady states - but does so, in contrast with the classical cases, with a speed proportional 
to Moreover, we show that close to the line of stable steady states, there is a 

correction to this large deviation approximation having a speed proportional to l/u. 
For a small but finite noise power u, P{A, B) is therefore the sum of two contributions: 
a large deviation approximation with speed 1 / which accurately describes the form 
of P{A, B) away from the attractor, and a large deviation approximation with speed 
1/z/, which describes P{A,B) close to the attractor. 

These competing large deviation terms are related in the AB model to different 
types of fluctuation paths having different large deviation speeds and, notably, to 
fluctuation paths of zero action, which we call "sub-instantons" . They are not related, 
we should mention, to the divergence of sub-leading prefactors entering in the large 



deviation approximation, discussed, e.g., by Berglund and Gentz 24,25 . We argue 



in the concluding section of the paper that these sub-instantons should arise in more 
general systems having multiple connected sets of steady states, and that they should 
lead, as in the AB model, to stationary probability densities having competing large 
deviation terms in the low-noise limit. 

The rest of the paper is organized as follow. In Sec. [2| we introduce the AB model, 
and then proceed in Sees. |3]j5]to obtain P{A, B) using three different approaches: a path 
integral approach which leads to results similar to those obtained in the framework of the 
Freidlin-Wentzell theory (Sec.|3]), a dynamical approach based on certain approximations 
of the AB model (Sec. |4]), and an approach based on the Hamilton- Jacobi equation 
(Sec. [5]). Throughout these sections, the analytical results obtained are compared with 
numerical results obtained by solving the Fokker-Planck equation directly. We conclude 
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Figure 1. Vector field of tfie AB model giving rise to a line ig of stable fixed points 
and a line of unstable fixed points. 

in Sec. [6] with some remarks on the generahty of our results. 
2. Model 

The model that we study is defined by the following set of two coupled (Ito) stochastic 
differential equations (SDEs): 

dA = {-AB - uA) dt + aA\fydWA 

dB = {J^ -vB)dt^OB^dWB, (3) 

where z/ > is a real coefficient that balances, together with the positive constants a ^ 
and Ob , the dissipation of the model and the intensity of the two uncorrelated Brownian 
motions WA{t) and WB{t). Henceforth, we refer to this model simply as the AB model. 
Note that the inclusion of the two dissipative or friction forces —uA and —uB prevents 
the system to escape (by diffusion) to 5 = oo as t — oo, and makes sure, therefore, 
that a stationary density P{A, B) exists. 

The zero-noise and zero-friction dynamics of the AB model is given by 

A = -AB 

B = A^. (4) 

This simple dynamics has the key properties that we referred to in the introduction, 
namely, 

(i) It has a continuous set of steady states, which corresponds here to the line A = 0; 

(ii) The dynamics is irreversible and converges to an attractor, corresponding here to 
the upper semi-line A = 0, B > 0, which we denote by ig- This semi-line is a line of 
stable steady states or stable fixed-points, as shown in Fig. [Tj The lower semi-line 
A = 0, B < 0, denoted by iu, is a line of unstable fixed-points||] 

I The connecting point (0, 0) is marginally stable. 
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Figure 2. Contour plots of P{A, B) obtained by direct numerical integration of the 
Fokker-Planck equation associated with the AB model for cta = = 1- The noise 
intensity v used in the integration is shown in the plots. Darker colors correspond 
to larger values of P{A,B). The region plotted for each noise power is different to 
highlight the concentration around the stable line. 

(iii) The AB dynamics conserves the quantity E = A? + B"^^ which we refer to as the 
energy. 

These properties are responsible for the competing large deviation scalings of 
P(y4, B) announced in the introduction. We shall show in the next sections that these 
scalings arise essentially from two very different fluctuation dynamics around Cg and 
which will be studied analytically. To support our results, we shall also present 
numerical results obtained by directly integrating the Fokker-Planck equation associated 
with Eq. ([3]) for different noise powers|§] Some of these results are shown in Fig. [2] as 
contour plots of B). We can already see from this figure that P{A, B) concentrates 
on the stable line £s as z/ — > 0, and that it is otherwise relatively isotropic away from 
is- Our analytical study of the AB model and of its fluctuations around is and iu will 
explain these two properties, among others. The main result that we shall obtain for 
c"A = o'b = 1 is the large deviation result 

lim-v^lnP(x) = J(A,5), (5) 

with 

I{AB) = ^-^{A' + B'f/\ (6) 

away from the stable line is- We shall also describe corrections of this large deviation 
result for finite values of u near is- 

§ The numerical integration was done with the routine NDSolve of Mathematica using double-digit 
accuracy and standard vanishing boundary conditions for square domains of the A-B plane. The 
domains used were [—5,5]^ for 1/ = 0.5 and v = 0.1, [—4,4]^ for — 0.05, and [—2,2]^ for ly = 0.025. 
Only a portion of these regions is shown in Fig. [21 
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3. Path integral solution 

We give in tliis section a first derivation of tlie large deviation form of P{A, B) following 



the classical theory of Freidlin and Wentzell 23 . This calculation is valid near the line 



Is of stable points and yields the large deviation speed Cy = 
3.1. Classical theory 

Consider a diffusion process {X{t)} in MP which is the solution of the SDE, 

dX(t) = /(X(t))dt + v^dW^(t), (7) 

where dW{t) are increments of the Brownian motion. It is known from the work of 
Freidlin and Wentzell (F-W) that the stationary density P{x) of this process satisfies, 
under certain conditions, a large deviation form or large deviation principle (LDP) in 
the low- noise limit — )■ oo, which we informally write as 

P(x) ^ e-^(")/", (8) 

to mean 

\mv-v\nP{x) = V{x). (9) 

I/— >-0 

Moreover, from the F-W theory, it is known that the rate function or pseudo-potential 
V{x) defined by this limit can be obtained from the following minimization problem: 

Vix) = inf inf L\x] (10) 

t>0 x{0)£O,x{t)=x 

which involves the Lagrangian or action 



L[x] = I C{x,x)ds, C{x , x) = -{x — f (x)) , (11) 







associated with a path {x(s)}*^o of process. The minimization in (10) is performed 
over all paths starting on the attractor O of the system at time t = and reaches the 
point X after a time t which is usually taken to go to infinity. 

The LDP for P{x) is akin to the semi-classical or WKB approximation of quantum 
mechanics, and can be explained heuristically as in quantum mechanics by expressing 
P{x) in path integral form: 

/•x{t)=x 

P{x) = lim P{x,t\x G 0,0) = lim / e"^''^]/'', (12) 

and by arguing that the probability to reach a point x is given, in the low-noise limit 
1/ — )■ 0, by the most probable path, called the optimal path or instanton, which starts on 
the attractor and reaches that point after a very long time. As this optimal path must 
have a minimal action under the terminal constraint x{0) G O and x{t) = x, we recover 



the result of (10). 



Freidlin and Wentzell 23 showed that this heuristic argument based on path 



integrals is rigorously valid for the SDE ([t]) provided essentially that O is a unique point 
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attractor of the deterministic dynamics x = f{x). If this dynamics admits many point 
attractors, then the result of Eq. (10) holds locally in regions G C that enclose 
single attractors and do not include characteristic boundaries, such as separatrices. 
Graham 17,26 also showed semi-heuristically that the F-W theory can be applied for 
non-point attractors, e.g., limit cycles or strange attractors. 

Here we apply the F-W theory to the AB model in a local sense by considering a 
region of the upper-half plane -B > surrounding the attracting line ig (see gray region 
in Fig. [T]). As we take the limit — )■ 0, we neglect the dissipative terms —uA and —uB in 
the full action of the model. Therefore, the action that we consider in the minimization 
problem is 

L[A,B]=fcds, C=-^{A + ABf + -^{B-A'f. (13) 
Jo 2(T^ 2a^ 

The reason for neglecting the dissipative terms in the action is that they are sub- 
dominant in the low- noise limit compared with the other terms, and lead to correction 
terms of order u in the rate function I {A, B). The net effect of neglecting these terms 
when solving the minimization problem of (10) is to look for instantons that emanate 
anywhere from the stable line ig of fixed points, which is the attractor of noiseless AB 
dynamics without the dissipation forces, i.e., Eq. (|4]), rather than emanating from the 
origin (0,0), which is the attractor of the noiseless AB dynamics with the dissipative 
forces. 

To find the instantons emanating from ig, we can solve the Euler-Lagrange 
equations 

d dC dC 

ds dA OA 
d dC dC 

where C is the Lagrangian density defined in Eq. (13), with the boundary conditions 
A{0),B(0) E is and A{t) = A, B{t) = B. By finding the solution of this equation for 
an increasing "hitting" time t, we then find a succession of approximations of the rate 
function I {A, B), which converge to I [A, 5) as t — )■ oo. 

For the AB model, the Euler-Lagrange equations read explicitly 

-\{A + AB- AB^) -\{A^- AB) = 



B 



-\{A^B + AA)-\{B-2AA) =0. (15) 
(^B 

This is a set of nonlinear, second-order equations of motion, which can be solved only 
numerically. In practice, it is easier to solve these equations by transposing them into 
an equivalent set of first-order Hamiltonian equations of motion defined by 

■ dH ■ on . m . on 

where 
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Figure 3. Red lines: Instantons reaching certain points {A, B) (black dots) from the 
stable attracting line of the AB model. Blue lines: Deterministic paths connecting the 
same point back to the stable attractor according to the noiseless dynamics of the AB 
model. Parameters: a a = <^b = 1- 



are the conjugate momenta associated with A and -B, respectively, and 

H = paA + pbB - C (18) 
is the Hamiltonian density. For the AB model, the Hamiltonian equations are explicitly 

A =a\pA-AB 
B =alpB + A^ 
Pa = PaB - 2pbA 

pB = PaA. (19) 

To numerically solve these equations, we use the fact that 7i is conserved in time 
and that H = on the attractor ig. As a result, instantons are such that H = for 
all times. Moreover, since the dynamics from the attractor is unstable, we numerically 
integrate the Hamiltonian equations backward in time instead of forward in time. This 
means that we start the integration from the point {A,B), and find initial conditions 
for Pa and pb such that the time-reverse dynamics of (19) leads to the attractor ig. 

A number of instantons obtained from this procedure are shown in Fig. [3] for 
different points (A, B) in the first quadrant of the A-B planejjj] The action associated 



with these paths, which yields the rate function I {A, B) according to Eq. (10), is shown 
in Fig. |4] as data points. The speed associated with this rate function is l/u. In the 
same figure, we show with the full and dashed lines the results of analytical calculations 
presented in Sec. |5} The agreement between the two sets of results will be discussed in 
more detail in these sections. 

An interesting property of the instantons, seen from Fig. [3} is that they are different 
from the natural decay paths of the system, i.e., the paths of the deterministic dynamics 

II By symmetry of the AB model, instantons for the second quadrant, i.e., on the left of must be 
the mirror images of the instantons of the first quadrant. 
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Figure 4. Rate function I{A, B) near is as a function of A for different values of B. 
Data points: Rate function obtained from the instantons shown in Fig. [Sj Dashed and 
full lines: Second- and fourth-order approximations, respectively, of /(A, B) obtained 
from the Hamilton-Jacobi method; see Sec.|5] Parameters: (ja = ctb = 1- 



starting from a point (A, B) away from the attractor and ending on the attractor (shown 
in blue). This difference is characteristic of SDEs that are not gradient, i.e., which cannot 
be written in the form 

X = -VU{x) + ^iit) (20) 

for a scalar function U{x). It is known that if such a form exists, then the rate function 
I{x) associated with the stationary distribution P[x) is simply given by I{x) = 2U{x). 
In this case, we moreover have that the instantons are the time-reversed version of the 
natural decay paths. 

For the AB model, these results do not apply: the SDE ^ is not of gradient- 
type, which means that we must explicitly solve the Euler-Lagrange equations or the 
equivalent Hamiltonian equations to find the instantons and their associated action. 
Physically, this also means that the AB model is a nonequilibrium system that violates 
detailed balance and gives rise to a non-vanishing probability current J = {J a, JbY i 
whose components are here given by 

va\ dP{A, B) 



J A = i-AB - i^A)P{A,B) 



2 OA 



J, = iA^-,B)PiA.B)-'^?^. (21) 

This current is actually a current loop, as shown in Fig. [5] it has the shape of a 
"squashed" double-solenoid and circulates through the origin. 



3.2. Sub -instantons 

The instantons shown in Fig. [3] are not strictly speaking the only possible instantons 
that the AB model admits. Interestingly, it is also possible to reach any points {A, B) 
from the attractor by following a fluctuation path constructed as follows (see Fig. [6]): 
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Figure 5. Streamlines of the probability current J for v — 0.05 and a a — = 1, 
with level colors representing The darker colors correspond to larger current 

magnitudes. 



(i) Start from any point on the stable attractor and go down to the origin; 

(ii) Follow the line iu of unstable fixed points down to the point (0,5') such that 

(iii) Follow the natural orbit of energy E = ^B'^ connecting (0, B') to (A, B). 

As each of these parts lies either on the line of stable and unstable points or follows a 
natural path of the deterministic system, the complete path thus constructed must have 
a zero action on the large deviation scale 1/z/. 

The existence of these zero-action instantons, or sub-instantons as we shall call 
them, gives us an indication that the dominant behavior of P{A, B) does not have a 
speed 1/z/; in fact, it has a speed l/y/i^, as we shall see below. The crucial point to note, 
however, is that P{A,B) has a sub-dominant LDP with speed 1/z/ near the attractor 
is (see gray region in Fig. [T]), which is precisely the LDP that we have obtained before 
with the F-W calculation. This sub-dominant LDP arises from instantons having a non- 
zero action emanating from is, whereas the dominant LDP, which scales with 1/a/z/, is 
related to fluctuation paths that are sub-instantons. The fact that these sub-instantons 
follow in the end the deterministic dynamics explains the observed isotropy of P{A, B) 
away from is, i.e., why P{A,B) has constant level curves along the vector field of the 
AB noiseless dynamics. The same indirect fluctuation paths also explain the existence 
of the probability current mentioned before. 

We show in the next two sections how to go beyond the F-W calculation to obtain 
the dominant LDP of P{A,B). Two methods are presented: the first is based on the 
solution of a time-scale separation of the solution of the AB model, whereas the second 
is based on a direct approximation of the solution of the Fokker-Planck equation. 
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Figure 6. Optimal paths. Tlie point in black in the A-B plane can be reached either 
by a direct, excited path (in red) with positive action or by an indirect path (in purple) 
that first goes along the stable and unstable lines of fixed points, and then follows a 
natural trajectory of the deterministic dynamics. The latter path has a null action. 



4. Dynamical analysis 

We analytically derive in this section LDPs for P(v4, i?) near 1^ and then near 1^ using 
the insight gained in the previous section. In treating and lu separately, we shall see 
how the stable and unstable dynamics, respectively, of each of these lines gives rise to 
different LDPs. As a prelude to these calculations, we obtain the stationary probability 
density of the energy defined by = + 5^ to show that it does not scale with i/, as 
is the case for weakly-perturbed Hamiltonian systems. 



^.i. Energy 

Given the original system of equations for A and S, displayed in Eq. (jsj), is found to 
evolve according to the Ito SDE: 

dE = -2z/Edt + z/(a^ + (T|)dt + 2v^(aAAdiyA + (TB5diyB). (22) 

This equation is not closed - it explicitly depends on A and B - but the equation for 
the mean energy is closed: 

i^ = -2(E) + (4 + a|). (23) 
From this, we find the stationary mean energy: 

(E), = (24) 

To find P{E), we solve the Fokker-Planck equation associated with the Langevin 
equation for E. For simplicity, consider first the case aA = ctb = <^ and let be any 
test function. Applying Ito's formula and averaging over the noises, we obtain 

i_mpi = <(,2 _ ^) + ^2 _ (25) 
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The associated Fokker-Planck equation for the energy pdf P{E, t) is thus 

2u dt dE > dE^ \ 

and admits the following stationary solution P{E): 

1 



P{E) 



(26) 



(27) 



It is easy to verify from this result that {E) = cr^, in accordance with the exact result 



of Eq. (24). Note also that P{E) does not depend on i/, as announced, which means 



that there is no LDP for P{E). 

For the case 7^ as, Ito's formula leads to 
Id (0(E)) /fal + al 



- E <P'{E) ) + {{alA' + alB') 0"(E)> , (28) 



2u dt \ V 2 

which is exact for any test function 0. By contrast with Eq. (|25|), the equation above is 



not a closed equation for E as the variables A and B are involved. However, because the 
typical values of the variable A are of order y^cii / \/B, which is much smaller than both 
\/E and B as soon as S> i^'^cf, it is natural in Eq. (28) to make the approximations 
B^ E and A"^ ^ E. One then obtains a closed equation for (j){E) whose associated 
Fokker-Planck equation is 



E 



a 



B 



1 dP _ d 
2^~dt ~ dE 
The stationary solution of this equation is 



P 



a 



(29) 



(30) 



where C is a normalization constant. We can verify that this expression reduces to 
the one we found in Eq. (27) when = o"b = o". It is also easily verified that the 



average energy obtained from this distribution is given by the exact relation displayed 



in Eq. (24). 



4-2. Dynamics around stable fixed points 

We now derive a local LDP for P{A, B) around ^s (see gray region in Fig. [T| using the 
fact that there exists a separation between the time with which A relaxes to (time of 
order 1/ B) and the time with which E varies (time of order l/v). Since E ^ B^ for A 
close to 0, the slow dynamics of E translates into a slow dynamics of B, which can be 
treated adiabatically with respect to A. Thus, we consider B to be constant and study 
the dynamics of the rapid variable A given by the first equation of the AB model: 

dA = {-AB - uA) dt + dWA (31) 

For fixed B, this equation is linear in A, which means that A evolves according to a 
Ornstein-Uhlenbeck process. In the stationary state, A is of order and the term vA 
can be considered negligible, and so we find 

BA^' 



P{A\B) 



B 



exp 



(32) 
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for the stationary distribution of the fast variable A given that B is fixed. This stationary 
distribution is reached on a time-scale of order 1/B. 
Next, we study the evolution of B given by 

dB = - B)dt + y^aB dWn, (33) 

using A' = Aj As we did for A, we can approximate A! to be an Ornstein-Uhlenbeck 
process with variance (A'^) = a\l{2B\ Accordingly, the fluctuations of the variable 
z/y4'^ must be such that {{uA''^ — (z/A'^))^) is of order z/^. Such fluctuations are much 
smaller than the effect of the white noise for B, and so, to leading order in u, the 
fluctuations of A'"^ can be neglected. The approximate dynamics for B to leading order 
in u is thus 

dB = iy(^^- B^dt + ^/J}aBdWB. (34) 

The associated Fokker-Planck equation is 

IdP _ d 
u~dt ~ dB 
and has for stationary solution 

P(5) = Cfi"^/"^e-^'/-|^ (36) 

where C is a normalization constant. By combining this solution with the solution found 
for P{A\B), we thus find 

P{A, B) = P{A\B)P{B) = C5^exp ^ ^ , (37) 

where C is again a normalization constant. The LDP extracted from this result has the 
form 

P(A, B) ^ e-^(^■^)/^ /(A, B) = ^, (38) 



2B 2 dB 



(35) 



a 



A 



and is nothing but the LDP associated with P{A\B). This shows that the large 
deviations of the AB model around the attractor are mainly the result of the diffusive 
fluctuations of A for B constant. 



The rate function /(A, B) of Eq. (38) is shown in Fig. |4|as a function of A and for 
different values of B as dashed lines. We see that there is a favorable comparison between 
this rate function and the rate function obtained from the F-W theory (represented again 
by the data-points). The difference with the F-W results can be explained by noting that 
the rate function derived in this section is only a second-order approximation of the rate 
function obtained by approximating the instantons as straight, horizontal paths joining 
a point (0, B) on to a point (A, B) near ^s. The correct instantons, as seen from 
Fig. [3} have a certain inclination or curvature in the B direction, which has the effect 
of reducing slightly the action of straight, constant-P paths. We shall see in Sec. [5] that 
the Hamiltonian-Jacobi method is able to capture this effect as a fourth-order correction 
to I{A,B). 
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4-3. Dynamics around unstable fixed points 

To obtain P{A, B) around the line of unstable fixed points, we need to construct a 
different approximation of the AB dynamics that accounts for the fact that trajectories 
starting at the origin have a small probability of diffusing down the line Cu and that, as 
soon as they depart from this line, they are rapidly expelled along the force lines of the 
AB dynamics, to reach the stable line Is'-, see Figs. [T] and [6j This dynamics creates a 
probability current because, eventually, the trajectories return to the origin by diffusing 
down the line Eg and start new loops by diffusing again down 

In order to qualitatively understand this fluctuation dynamics, we consider the 
following approximation of the AB model: 

dA = aAAdt + aAv^dWA 

dB = -uB dt + asV^dWe, (39) 

which involves a "top-hill" diffusion process for A (a^ > 0) and a "down-hill" diffusion 
for B. The equation for B is the exact dynamics of B on the line A = 0, whereas the 
"top-hill" term oaA is analogous to the term —AB in the original AB model. For now, 
we consider oa to be independent of B in order to have a solvable model which will 
help us gain a qualitative understanding of the unstable dynamics. We discuss the case 
aA = —B at the end of this section. 



For aA constant, Eq. (39) describes an unstable linear Gaussian process. To find 
a stationary distribution P{A, B) for this process, we have to consider that there is a 
probability (or current) source at the origin O. In this case, the probability to reach a 
point {A, B) in the lower plane is given by integrating the probability P(A, B, t\0) that 
a trajectory starting at the origin O reaches the point (A, B) after a variable time t, to 
which we multiply the stationary probability density of starting at O: 

POO 

P{A,B) = P{0) PiA,B,t\0)dt. (40) 
Jo 



For the decoupled dynamics of Eq. ( p39j ), the propagator P{A, B,t\0) is simply given by 

(41) 



P(A,B,t\0) = exp 

27i^{A^){B^) 



A' B' 
2(A2) + 2(^2) 



where 



and 



(^2^ = ^(e2"At_i) (42) 
2aA 

2 

(^2) = ^(l-e-2^"^*). (43) 

The difference between the behavior of the variance of A and B plays an important 
role in determing the LDP of P{A, B) in the low-noise limit. The variance of A grows 
exponentially for large times, as a result of the instability of lu-, whereas that of B goes 
to a constant on a long time-scale of order t ~ {2a\v)^^. Taking the limit — )■ 0, and 
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anticipating that the integral will be dominated by contributions with times such that 
ut <^ 1, we obtain 

POO 

P{A = 0,B)^C / e-'-At-By{2alut) ^44^) 
Jo 

for the probability density on with C a normalization constant. This can be put in 
the form of a Laplace integral with the change of variables t' = a\vt\ 

POO 

P{A = 0,B)^C / e-5(*')/v^ dt' (45) 
Jo 

where 

In the limit z/ — )■ 0, the integral is therefore dominated by its maximum integrand, which 
is here located at t' = B/y/2aA- As a result, we find 

P{A = 0,B)^ e-'^(^)/^, ifi{B) = ^^?^B. (47) 

Note the change of speed in this LDP: l/v^ instead of l/u. 

From this simple "top-hill" linear model, we conclude that the low probability 
P(0, B) for the system to be on the unstable line iu results from the balance between 
the probability to reach a finite value B after a time t and the probability that the system 
is not expelled by the unstable dynamics before this time t. The rare coincidence of these 
two conditions gives the large deviation result. The preceding computation suggests that 
the leading contribution corresponding to those rare events correspond to an optimal 
time t{B) proportional to y/u. 

This qualitative behavior also correctly describes the exact AB dynamics with 



aA = B. In this case, one can explicitly solve the two equations shown in (39) with 
aA = B. However, we have found more convenient to use in this case an asymptotic 
expansion of the Hamilton- Jacobi equation, which we discuss in the next section. For 



now, we just note that substituting ua = B 'm. Eq. (47) leads to ^{B) oc B^/'^ /aB-, a 
result that will be confirmed in the next section. 

To close this section, note that the density P{A, B) away from is obtained by 
radially radiating the value of P{0,B) obtained above, i.e., by setting P{A',B') = 
P(0, B) for all {A', B') such that A'"^ + B'"^ = B"^. This is expected from the numerical 
results shown in Fig.[2]and follows, more precisely, by noticing that trajectories departing 
infinitesimally from are radiated radially by the vector field of the AB model. As this 
instability results only from the deterministic part of the AB equations, the probability 
density must have circular level curves. 

5. Hamilton- Jacobi approach 

We give in this section a different derivation of the LDP of P{A, B) near ^s and 
obtained by solving the stationary Fokker-Planck equation of the AB model, which 
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reads 

A(ABP) - A(^.p) ^ , ,|_(BP) ^ ^ . 0. (48) 

As before, we look for solutions having the LDP form 

P{A, B) = (f){A, B) e-^-^^^'^), (49) 

where I{A, B) is the rate function and 0(A, B) a function of A and B that does not 
depend on u. Moreover, following the previous section, we choose the noise scaling 
Ci, = l/u around Is and Cy = 1/v^ around £u- 

The advantage of the Fokker-Planck method is that it yields more precise 
expressions for the rate function in the form of series expansions, which are not based, 
as in the previous section, on approximations of the AB model. However, to be able to 
use this method, we need to know the correct large deviation speed of P{A, B), i.e., the 
correct scaling of c^, with u. 



5.1. Stable fixed points 



Inserting the ansatz (49) with Cy = l/p m. the Fokker-Planck equation (48) yields, to 



order l/i^, the following partial differential equation: 

This is a Hamilton- Jacobi equation corresponding to the weak-noise limit of our model. 
Unfortunately, we are unable to solve this equation explicitly. However, we expect the 
solution to be localized close to A = 0, and so it is natural to look for an expansion of 
I for small A. By symmetry of the AB model, the dependence of P{A, B) on A must 
be even. Thus, assuming an expansion of I of the form 

liA, B) = Io{B) + h{B)A^ + h{B)A' + 0{A^), (51) 



in Eq. (50), we directly obtain 



/(A B) = ^A^ - + 0{A% (52) 



Figure |4] shows how this result compares with the calculation of the rate function 
based on the instanton approximation presented in Sec. |3} The full colored lines on this 



plot show I (A, B), as given in Eq. (52), as a function of A for different values of B. The 
dashed colored lines represent the same result but truncated to second order. The data 
points, as mentioned in Sec. [3} are the results of the instanton approximation. The near 
perfect match between the data points and the full lines confirm the consistency of the 
Hamilton- Jacobi and instanton results. 

The result of Eq. (52 ) is also compared in Fig. [Tjwith results obtained by numerically 
integrating the Fokker-Planck equation. We see in this figure that, for small values of 
the numerical results are relatively close to /, at least for points [A, B) near ig. It must 
be said that it is rather difficult to obtain reliable results by numerically integrating the 
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Figure 7. Numerical test for I{A,B). Black dashed line: Second-order result for 
I{A, B) obtained from the Hamilton-Jacobi method. Black full line: Fourth-order 
analytical result for /(A, B) obtained with the same method. Colored lines: Rate 
functions obtained from the numerical integration of the Fokker-Planck equation of 
the AB model for different values of v. Parameters: ct^i = ctb = 1. 



Fokker-Planck with small noise intensities, especially away from the attractor. In our 
case, we could not obtain reliable results for values of v smaller than about v = 0.025. 
For this reason, the numerical results of Fig. [7] are presented to check that our analytical 
results have the correct scaling in z/, and that the rate functions obtained numerically 
and analytically are qualitatively similar. 

To obtain the correction factor 0, we need to consider the Hamilton-Jacobi equation 
for (f) obtained at order of the Fokker-Planck equation: 

b^ + ab£-A'^-a-4.-b-4. 

"^^"''" + 2^-^(1^^ +2^=0. (53) 



2 V^^^ dAdAj 2 V^^^ OBdB^ 
As before, we attempt to solve this equation by expanding in even powers of A. Thus 
we assume 

B) = MB) + MB)A^ + O(A^), (54) 

and find at order A"^, 

K<^\ + <^l)B(t)'o + (12^2 - Qa\ - 3(t|)0o + S^V^ 0i = 0, (55) 

where 4>q{B) is the i?-derivative of (f)o{B). This equation is not a closed differential 
equation for 0o as it involves (pi. To get a closed set of equations, we find a 
complementary equation by looking at the Hamilton-Jacobi equation obtained at order 
i/, which reads 

<t\ ( d^I ^ dl dip d\ 



r^ + 2 



2 V^^^ dAdA dA^ 
al f d^I dl dip d\ , 
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Figure 8. Numerical test for I{A,B) on tlie line B = —1. cr^ = o"b = 1- Black 
line: Analytical result for I {A, B). Colored lines: Rate function /(A, B) obtained from 
the numerical integration of the Fokker-Planck equation of the AB model for different 
values of v. 



We seek a solution of the form ip{A,B) = (po{B) + 0(^4^) for the function carrying 
the i/-order correction. Substituting in the equation above the expressions of / and (j) 



shown in (52) and (54), respectively, we obtain at order A^: 
al^l + 250'o + 400 + 2a\(t)i = 0. 



(57) 



With this equation and Eq. (55), it is possible to eliminate the function (pi to obtain a 
closed differential equation for 0o: 



4SVi,0'o'+[85^-45(al + a^) 



6a\ 



3a 



Bm 



0. 



Looking for solutions of the form CB"e we then find 



CB^e-^ /''B, 



(5J 



(59) 



where C is a normalization constant. This is the only normalizable solution to Eq. (58) 



I.e.. 



With the solution of / truncated to order 0{A'^) and truncated to order 0{A^), 

\2\ 



PiA,B) 



')o + 0(A ), we can finally write down our approximation for P{A, B): 
B^' 



B^B 



exp 



X exp 



BA"' 



B 



+ 



+ 0{A') 



(60) 



which reproduces the result of Eq. (37) to second order in A in the exponential. 



5.2. Unstable fixed points 

Knowing from the previous sections that P{A, B) is isotropic around the line of 
unstable points, it makes sense to solve the Fokker-Planck equation in polar coordinates: 

dP d ^ , 2 ■ 2 2 2 

— z/r+— (cT^sm y + o"^ cos t 



dt 



dr . 



2r 
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+ 
+ 

+ 



9_ 

2dr^ 
V 92 



r 008 6* 



p 



[(a^cos2^ + (T|sin2^)P] 



2drde 
2W 



sin 2^ 



-ex, 



sin 20 



P 



sm 



2 cos^ Q 



P 



Ga = CTb 

equation 



(61) 

At order z/~^/^, and setting 
1 for simplicity, this equation reduces to the simple partial differential 



using the large deviation ansatz P{A, B) 



— r cos 9 



m 

39 



0, 



(62) 



which shows that / does not depend on 9 for A ^ (i.e., 9 7^ ±7r/2). To find the 
r-dependence of /, we have to consider the Hamilton- Jacobi equation associated with 
the 0th order of the Fokker-Planck equation in z/, which reads 



d_ 
89 



r 0cos( 



+ 



1 /dJ\ 



2 V dr / 



0. 



(63) 



We look for a solution of this equation close to 9 = —njl by making the change of 
variables 9 = — 7r/2 + u and by defining 0(m, r) = 0(— 7r/2 + u, r). The equation above 
then becomes 

2 



d_ 
du 



r sm u 



dl 

dr 



0. 



(64) 



By symmetry, (p must be even in u, which implies that dcp/du = at u = 0, and so we 
find 

2" 

'0(O,r)=O. (65) 



dl 

dr 



Assuming that 0(0, r) > 0, we then obtain 



/(r) 



2^2 



.3/2 



as the solution for Jj^ In terms of A and B, this yields the LDP 



P{A,B) ^ e 



-/(A,iJ)/v^ 



where 



I{A,B) 



2V2 



{A^ + B'f\ 



(66) 
(67) 
(68) 



This rate function is different from the one obtained in Sec. HI This is because the 
dynamical model studied in Sec. |4] is an approximation of the AB model in which the 
dynamics of A is decoupled from B. This approximation is not assumed here, so we 

% Note that / must be positive and must vanish on the stable fixed point (0, 0). 
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Figure 9. (a) Fluctuation paths for disconnected sets of stable (As) and unstable 
(Au) fixed points. In this case, any path reaching a point x must contain an instanton 
part, even if it goes through A^ (blue path), so its large deviation probability must 
have a speed l/i^. (b) Fluctuation paths for connected sets of stable and unstable fixed 
points. In this case, there are fluctuation paths that reach x via "sub-instantons" (blue 
path) instead of instantons (black path). Sub-instantons have a large deviation speed 
smaller than 

should take the solution above as a more precise expression of the rate function. The 
essential point to notice is that both approaches confirm the speed l/v^ of the LDP of 
P{A, B) near which also appears to qualitatively match the scaling of the numerical 
Fokker-Planck results; see Fig. |8j As before, it is difficult to push the numerical scaling 
analysis beyond u = 0.025, as the numerical integration of the Fokker-Planck equation 
becomes unstable for smaller values of u. 



6. Discussion 

The two LDPs that we have derived for P{A, B) can be thought of as an expansion of 
P(y4, B) that keeps its two first dominant terms, so that 

P{A, B) = Ci e-^(^'^)/^ + C2 e--^(^'^)/^ (69) 

As discussed, the dominant term in this expansion is the \/v term associated with the 
repelling line and is "non-classical" in the sense that it does not arise from instanton 



fiuctuation paths that are characteristic of the theory of Freidlin and Wentzell 23 . The 
second term with large deviation speed 1/z/ is the term that arises from this theory: it is 
associated with instantons emanating from the attracting line and becomes important, 
for small but finite in the vicinity of this line. 

The presence of the stable and unstable lines of fixed points is important for 
obtaining these results, but what is also important is that these two sets of fixed points 
are connected and that the system can diffuse from one to the other. If the attracting 
and repelling lines were not connected, then all the fluctuation paths reaching a given 
point away from these lines would necessarily be instantons having a non-zero action 
at the 1/z/ order in the exponential of P(v4, B\ which means that P(v4, i?) would then 
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globally satisfy an LDP with speed l/u. This is illustrated in a general way in Fig. [9] 
There we see that a point x can be reached from an attracting set Ag either directly 
with a single instanton (black path in Fig. [9^), or with a fluctuation path that reaches 
the repelling set Au before reaching x (blue path). The second part of latter fluctuation 
path, going from A^ to x, follows the natural dynamics of the system, and so carries zero 
action, while the first part, going from Ag to Au, is an instanton with non-zero action 
at the scale l/u. Consequently, the whole path must inherit the action of the instanton, 
which leads us to conclude that large deviations in this scenario are instanton-related 
and, thus, have the speed l/u. 

If the sets Ag and A^ of fixed points are connected, as is the case for the AB model, 
then the part of the fluctuation path that goes from the attracting set to the repelling set 
(see Fig. |9}d) is not an instanton anymore, at least not an instanton with positive action 
at the scale l/u. In the AB model, it is a simple diffusing path or suh-instanton that 
contributes, as we have seen, to a factor in the stationary probability density, 

where c is some positive constant. In general, this is greater than the probability of 
"true" instantons, which scales like e"'^/'' as — t- 0. It might be the case that, for points 
close to the attracting set, and for small but finite values of z/, instantons can have a 
probability greater than the probability of sub-instantons. This happens, for example, 
in the AB model near the attracting line. However, the fact is that, because of their 
lower speed, sub-instantons necessarily have a probability higher than the probability 
of instantons in the limit z/ — )■ 0, which means that they determine the dominant LDP. 

We expect this scenario to apply to more realistic systems having extended and 
connected sets of steady states. In general, the probability of diffusing paths or sub- 
instantons going from an attracting set to a repelling set is likely to depend on the 
particular system considered, but what should be clear is that these paths have, in the 
low- noise limit, a probability greater than the probability of instantons. As a result, 
they must determine, if they exist, the dominant large deviation form of the stationary 
probability density of the system, when such a stationary density exists. In the future, 
it would be interesting to determine minimal conditions for the appearance of sub- 
instantons or related fluctuation paths, and to see whether these generally arise in 
stochastic systems violating detailed balance, as seems to be suggested by our results. 
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